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The  spin-flop  transition  in  a  uniaxial  antiferromagnet  defines 
three  critical  fields;  that  of  the  true  thermodynamic  transition  and 
those  limiting  the  local  stability  of  antiferroraagnetic  and  spin- 
flop  phases.  The  latter  two  are  calculated  by  a  spin-wave  analysis. 
The  results  are  in  qualitative  agreement  with  the  hysteresis  observed 
by  Schelleng  and  Friedberg  in  the  spin-flop  transition  in  MnBr2.4HgO. 
The  spin-wave  spectrum  of  the  spin-flop  phase  is  given,  and  the 
magnetization  in  the  spin-flop  phase  is  found  to  increase  more  rapidly 
than  linearly  with  field  because  of  quantum  corrections  absent  in 
molecular  field  theory. 
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We  consider  a  uniaxial  antiferromagnet,  at  or  very  near  zero 
temperature,  with  an  applied  field  parallel  to  the  unique  easy  axis. 
As  the  field  is  increased  a  transition  occurs  to  a  spin-flop  phase, 
in  which  the  spins  are  almost  orthogonal  to  the  field.  Such  tran¬ 
sitions  have  been  observed  by  many  investigators,  and  in  certain 
materials  considerable  hysteresis  may  be  observed  in  the  transition, 

with  different  critical  fields  as  the  transition  is  traversed  up- 
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ward  and  downward  in  field. '  '  '  This  is  to  be  expected  on  the 

following  basis.  As  the  field  is  increased  to  a  critical  value 

Hu,  a  local  instability  occurs  when  one  of  the  spin  wave  frequencies 
c 

of  the  antiferromagnetic  configuration  becomes  negative.  As  the 

field  is  decreased,  a  local  instability  in  the  spin-flop  phase 

4 

occurs  at  a  lower  field  H  ,  when  one  of  the  spin  wave  frequencies 

C 

of  the  spin-flop  phase  becomes  negative.  Between  the  two  critical 

fields  and  H*,  there  exists  a  third  field  H°,  at  which  the  free 
c  c  c 

energies  of  the  two  phases  are  equal.  Whether  the  actual  observed 

transitions  occur  at  H°,  without  hysteresis,  or  at  and  H* , 

c  c  c 

with  hysteresis,  depends  upon  the  presence  or  absence  of  nucleating 
centers  and  local  inhomogeneltles.  The  problem  is  fully  analogous 
to  the  super-heating  and  super-cooling  metastability  in  a  con¬ 
ventional  gas-liquid  first-order  phase  transition. 

An  alternative  mechanism  for  the  hysteresis  has  been  suggested 
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by  Date  and  Hagata.'  '  They  ascribe  the  hysteresis  to  the  exist¬ 
ence  of  a  fourth-order  anisotropy,  in  addition  to  the  dominant 
uniaxial  anisotropy.  Tills  fourth-order  anisotropy  introduces 


energy  minima  along  both  the  longitudinal  and  transverse  axes,  with 
an  intermediate  energy  maximum.  The  hysteresis  is  associated  with 
the  difficulty  of  overcoming  this  intermediate  energy  barrier. 
However,  the  hysteresis  certainly  occurs,  and  is  a  large  effect, 
even  in  materials  with  very  small  fourth-order  anisotropy,  as  in 
MnBr2.4H20. ^  * 

For  simplicity,  we  assume  a  model  decomposable  into  two  equiva 
lent  sublattices,  with  the  nearest  neighbors  of  a  spin  on  one  sub¬ 
lattice  lying  entirely  on  the  other  sublattice.  The  anisotropy 

z  2 

is  taken  simply  as  -K(S  )  for  each  ion,  and  2S  is  the  exchange 
interaction  of  nearest  neighbor  pairs.  Labelling  the  spins  on 
one  sublattice  by  f  and  those  on  the  other  sublattice  by  g,  the 
Hamiltonian  is 


Spin  wave  analyses  of  the  antiferromagnetic  phase  have  been 
given  by  P.  W.  Anderson^  ^  and  by  R.  Kubo^,  to  which  we  make 

one  amendment.  In  application  of  the  Holstein-Primakoff  trans- 

z  2  +2 

formation  the  K(Sf)  operator  is  replaced  by  KlS-a^a^)  * 

K[S2-2Sa£af  +  (a*af)2]  where  a+af  is  a  boson  number  operator. 

The  spin  wave  theory  can  be  looked  upon  as  an  expansion  in  powers 

.  2 

of  the  spin  deviation,  and  (a^a^)  can  then  be  neglected.  Alter¬ 
natively  all  operators  can  be  replaced  by  simple  operators 
which  give  the  correct  matrix  elements 
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in  the  low-lying  states  (which  alone  are  important  in  the  spin  wave 
regime)  even  though  relatively  large  errors  are  introduced  in  the 
high-lying  states;  this  is  the  procedure  suggested  by  Van  Kranendonk 
and  Van  Vleck.^  Thus, for  the  lowest  two  levels(in  which 
afaf  "  0  and  1)  we  have  (afaf)2  *  afaf*  and  ln  Procedure 

K(S*)2  ■  K[S2  -  ( 2S-1 )  a^a^].  We  note  that  for  S  »  1/2  both  the 
right-hand  and  left-hand  members  then  reduce  to  K/4,  whereas  the 
equation  is  inconsistent  if  (a+af)2  is  neglected.  The  alteration 
in  the  latter  formalism  is  made  by  replacing  2SK  by  (23-1 )k,  or 
K  by  C2K,  where 


The  familiar  result  of  the  analysis  is  that  the  spin  wave  spectrum^'^*-^ 
of  the  antiferromagnetic  phase  is  composed  of  a  doubly  degenerate 
branch  with  an  energy  gap,  concave  upward  at  k  *  0,  and  approaching 
linearity  with  increasing  k.  The  applied  field  splits  the  degeneracy, 
driving  one  branch  upward  and  one  downward.  At  the  critical  field 
the  k  m  0  mode  of  the  lower  branch  is  driven  down  to  zero  fre¬ 
quency,  and  this  occurs  at  a  field  which,  with  the  i2  correction 
Inserted,  is 

where  z  is  the  number  of  nearest  neighbors. 

L 

To  calculate  Hc  we  formulate  a  spin  wave  theory  of  the  spin- 
flop  phase.  We  choose  coordinate  systems  Xj,  y,  z ^  and  x2,  y,  zg 
to  characterize  spins  on  the  f  and  g  sublattices  respectively. 
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The  z !  and  z2  axes  lie  in  the  first  and  second  quadrants  of  the 
x-z  plane  respectively,  and  each  makes  an  angle  6  with  the  z-axis. 

In  choosing  these  coordinate  systems  it  is  anticipated  that  the 
spins,  initially  in  the  x-z  plane,  make  only  small  oscillations  in 
the  vicinity  of  that  plane.  In  particular,  they  do  not  precess  around 
the  z-axis,  away  from  the  x-z  plane.  To  ensure  that  this  is  so  we 
add  an  additional  planar  anisotropy  to  the  Hamiltonian; 

W*  14  [£  o> 

At  the  end  of  the  analysis  we  shall  show  that  the  planar  aniso¬ 
tropy  has  no  essential  effect,  that  the  spectrum  and  critical 
field  are  very  weak  functions  of  K1,  and  that  we  can  take  »  0 
without  complication. 

We  make  the  Holstein-Primakof f  transformation,  followed  by 
the  replacement  by  operators  appropriate  to  the  low-lying  states 
as  discussed  above. 

=  S  -  al  a(  (5) 


(6) 


(7) 

and  similarly  for  S+,  s”.  These  operators  have  correct  matrix 

6  6  6 

elements  among  the  three  lowest-lying  states  of  each  spin.  Sub¬ 
stituting  in  the  Hamiltonian  H’ ,  we  neglect  terms  whose  matrix 
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elements  would  be  quadratic  in  the  amplitude  of  the  second  excited 
state,  such  as  a^a^a^a^,  or  those  whose  matrix  elements  would  be  pro¬ 
portional  to  the  product  of  the  amplitudes  of  the  first  and  second 
excited  states,  such  as  a^a^a^.  A  term  such  as  afafafafafaf  ls 
replaced  by  2a+af.  We  thereby  obtain  a  Hamiltonian  involving  zero, 
first,  and  second  order  terms  in  the  boson  operators.  The  condition 


that  the  first  order  terms  vanishes  determines  0: 

N 

^  X  2  1 

The  resultant  Hamiltonian  is  most  conveniently  written  in  terms  of 


spin  wave  operators,  defined  by 


and 


(10) 


whence 

»':E.  +  *£(<£«< +1*1^)  + 


(ID 


I 
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where 


A  *  ails 

=  ~2*XS  <1^4  ^ 
-  2*X5c«>**  ^ 


(12) 


(13) 

(1M 

(15) 


and  is  defined,  as  usual,  by  the  summation  of  exp  (ik.6^) 
the  z  nearest  neighbors  at  positions  6  The  Hamiltonian  is 
diagonalized  by  the  transformations 

S  -i-NvKWv**) 


where 


6*  * 


over 


(16) 

(17) 


x,  -  (L'Ml'*  •**)*' 

\  a  -eA  ) 

M-* B-J^^  <"* 
U«a+d4»eA) 


(18) 


(19) 


W'*F.  (20) 


giving 
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with  the  spin  wave  energies 


and  with  K  and  Kj  denoting  the  "reduced”  anisotropy  constants 


k«k.A»»t> 


It  will  also  be  recalled  that  5 


>  £,  *  *,  /M) 

] 1-(2S)-1 


(23) 


and 


2  -1 

cos  6  *  ^H(2S(2zJ  -5  K)]  ,  in  terms  of  which  the  spectrum  of  the 

spin-flop  phase  is  completely  defined. 

For  small  anisotropy  the  UJ1(k)  branch  is  almost  linear  in  k, 
because  of  the  first  square  root,  whereas  the  “>2(k)  branch  has  a 
large  k  =  0  intercept.  Consequently  we  refer  to  the  “^(k)  branch 
as  the  acoustical  branch  and  the  *2(k)  branch  as  the  optical  branch. 
The  optical  branch  has  vanishing  initial  slope,  a  negative  initial 
curvature  for  large  applied  fields  (cos  wfc  *  ***  K1  “ 

and  a  positive  initial  curvature  for  small  applied  fields,  as  in¬ 
dicated  in  Fig.  1.  The  change  in  Initial  curvature  occurs  because 
the  values  of  *2(k)  for  small  k  decrease  rapidly  with  increasing  H, 
whereas  the  decrease  is  less  rapid  for  large  k.  In  contrast  the 
applied  field  Increases  the  frequency  of  modes  in  the  acoustical 
branch. 

L 

As  the  field  is  decreased  to  the  critical  value  H  the  fre- 

c 

quency  of  the  k  =  0  optical  mode  vanishes  (see  Fig.  1).  This 
critical  field  is  found  by  taking  Y,  =  1  and  setting  the  second 
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radical  In  equation 

Co1  = 


(21)  equal  to  zero,  giving 

-  to-nfr 


or 


S'K) 


tow* -t(bt)rt 

Hit  +f(M-r)K 


(24) 


(25) 


We  now  note  that  5-1  for  large  spin,  and  in  the  classical 
limit  ( S*to)  the  critical  field  becomes  independent  of  K^.  For 
even  small  spin  the  fictitious  anisotropy  has  little  effect;  in 
the  most  extreme  case  ( S=»l )  the  coefficient  of  is  only  .21, 
whereas  the  coefficient  of  K  is  1.20.  The  effect  of  the  aniso¬ 
tropy  Kj  for  small  spin  values  is  traceable  to  the  quantum 
fluctuations  of  S? ,  which  enables  these  spins  to  sense  the  energy 
out  of  the  x-z  plane.  In  accordance  with  the  reality  of  the 

physical  situation  we  henceforth  take  Kj  =  0. 

ft 

For  MnF2  the  ratio  K=  K/(2zJ)  is  m  .01,  so  that  the  two 

ii  A 

fields  and  Hc  differ  by  only  about  2%.  But  for  MnBr2.4HgO 

the  field  2s(2zj)/U  is  approximately  25  kilo-oersteds,  whereas 

u  A 

2SKA1  is  approximately  5  k-oe.  Hence  Hc  and  Hc  differ  by  about 

2  k-oe.,  with  a  mean  of  about  8  k-oe.  This  in  general  agreement 

with  the  scale  of  the  hysteresis  observed  by  Schelleng  ar.d  Fried- 

berg^  in  the  spin-flop  transition. 

It  should  perhaps  be  mentioned  that  the  field  can  be  found 

c 

alternatively  by  extrapolating  to  zero  frequency  the  field  H(“) 
required  for  antiferromagnetic  resonance  at  frequency  w.  Thus  the 


11 


separation  of  the  critical  fields  and  H°  can  be  observed  in 
principle  even  in  materials  in  which  impurities  destroy  the  meta¬ 
stability  of  the  ''superheating''  or  'supercooling''  states. 

The  net  matnetization  ^  S^and  the  sublattice  magnetization 
^  are  of  interest.  These  quantities  are  i*elated  by 


x  **  y  c'*  \ 


(26) 


and  the  sublattice  magnetization  is 


<?>  -  s*t  -^l{— 


(27) 


Here  n^k)  is  the  average  occupation  of  the  acoustical  mode  of  wave 
vector  k;  n1(k)  *  £ exp  Oh*^ ( k ) ] -l^j  _1,  and  similarly  for  n2(k)* 
Several  limiting  cases  of  interest  can  be  identified  im¬ 
mediately.  Let  T  =  0  (or  p  =  “^  H  =  0,  and  K  =  0.  we  then  are 
discussing  the  ground  state  of  the  simple  antiferromagnet,  and 
equation  (27)  reduces  directly  to 

(sl">  *s  "t  £  “'I  < 

which  is  the  familiar  result  of  Anderson.  By  numerical  cal¬ 
culation  of  the  summation  in  equation  (28)  Anderson  has  found  a 
reduction  of  approximately  . 08  in  the  effective  spin  from  its 
'ideal'  value  of  S. 


In  the  opposite  limit  we  again  take  T  =  0  and  we  note  that 
when  H  takes  the  critical  value 

=2S(«I-C‘K) 

then  0=0  and  ^S^1)  =  S.  Thus  the  spin-flop  phase  undergoes  a 
second  order  phase  transition  to  a  fully-aligned  phase,  with 
fully  saturated  spin  components.  This  second-order  transition  has 
been  discussed  by  Palk^^  for  the  special  case  of  zero  anisotropy 
by  a  variational  method,  and  by  Anderson  and  Callen^  by  a  Green 
function  method  which  extends  the  theory  to  higher  temperatures. 
The  latter  authors  have  also  analyzed  the  critical  field  in 
the  intermediate  temperature  region,  and  the  thermodynamics  of  the 
antiferromagnetic  and  paramagnetic  phases. 

As  the  field  is  Increased  at  constant  temperature,  through 
the  spin-flop  range ,^S^lncreases  through  two  contributions. 
Firstly,  cos  0  increases  linearly  with  H,  and  secondly  the  sub- 
lattice  magnetization  \  S^, J  /  increases  (becoming  equal  to  S  only 
at  HQ).  Consequently  increases  faster  than  linearly  with  H, 

as  contrasted  with  the  strictly  linear  variation  predicted  by 

(8) 


molecular  field  theory 
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